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Abstract 

We show - by explicit computation of first-order corrections - that the QCD fac- 
torization approach previously apphed to hadronic two-body decays and to form 
factor ratios also ahows us to compute non-factorizable corrections to exclusive, 
radiative B meson decays in the heavy quark mass hmit. This removes a major 
part of the theoretical uncertainty in the region of small invariant mass of the 
photon. We discuss in particular the decays B — > K*^ and B — > K*£~^i~ and 
complete the calculation of corrections to the forward-backward asymmetry zero. 
The new correction shifts the asymmetry zero by 30%, but the result confirms our 
previous conclusion that the asymmetry zero provides a clean phenomenological 
determination of the Wilson coefficient Cg. 



(submitted to Nucl. Phys. B) 



1 Introduction 



Radiative B meson decays are interesting because they proceed entirely through loop 
effects. The chiral nature of weak interactions implies additional suppression factors 
thus enhancing the sensitivity of radiative decays to virtual effects in theories beyond the 
standard theory of flavour violation, if these new effects violate chirality. The observation 
of 6 — s> S7 transitions with a branching ratio of about 10~^ (in agreement with the 
standard theory) has since provided significant constraints on extensions of the standard 
model. The semi-leptonic decay b —>■ si~^i~ is expected to occur with a branching fraction 
of a few times 10~^ and will probably be first observed in the near future. 

Radiative flavour-changing neutral current transitions can be measured inclusively 
over the hadronic final state (containing strangeness in the particular case considered 
here) or exclusively by tagging a particular light hadron, typically a kaon. The inclusive 
measurement is experimentally more difficult but theoretically simpler to interpret, since 
the decay rate is well and systematically approximated by the decay of a free b quark 
into light quarks and gluons. In the long run, however, the easier detection of exclusive 
transitions requires the development of a systematic theoretical framework for these 
modes as well. 

The theoretical difficulty with exclusive decay modes is usually phrased as the need 
to know the hadronic form factors for the B — > K^*^ transition, but this is only one 
aspect of the problem. Even if the form factors were known with infinite precision, the 
present treatment of exclusive, radiative decays would be incomplete, because there ex- 
ist "non-factorizable" strong interaction effects that do not correspond to form factors. 
They arise from electromagnetic corrections to the matrix elements of purely hadronic 
operators in the weak effective Hamiltonian. We compute these non-factorizable correc- 
tions in this paper and demonstrate that exclusive, radiative decays can be treated in a 
similarly systematic manner as their inclusive counterparts. As a result we obtain the 
branching fractions for B — K*'~f and B — > K*i^i~ for small invariant mass of the lep- 
ton pair to next-to-leading logarithmic (NLL) order in renormalization-group-improved 
perturbation theory. (In another commonly used terminology that takes into account 
the l/os enhancement of the semileptonic operator Og the result would be referred to 
as next-to-next-to-leading logarithmic. However, in order not to have to distinguish in 
terminology the photonic and semileptonic decay, we shall refer to both as NLL when 
actually NNLL is meant for the semileptonic decay. Furthermore, since some of the 
3-loop anomalous dimensions required for the semileptonic decay are still not computed, 
the accuracy is not strictly NNLL.) 

Our method draws heavily on methods developed recently for other exclusive B de- 
cays. In Q power counting in the heavy quark mass limit and standard factorization 
arguments for hard strong interaction processes have been used to demonstrate that 
decay amplitudes for hadronic two-body decays of B mesons can be systematically com- 
puted in terms of form factors, light-cone distribution amplitudes and perturbative hard 
scattering kernels. A similar reasoning has been applied by two of us ^ to hadronic 
form factors for B decays into light mesons in the kinematic region of large recoil of 
the light meson. We also noted that a combination of the work on non-leptonic decays 
and on form factors could be used to compute non-factorizable corrections to radiative 
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decay amplitudes. Furthermore, the ten different form factors that exclusive decays may 
depend on can be reduced to only three p. More precisely (but still schematically), the 
amplitude can be represented as 

{trKl\H,^\B) = Caia + <fB ® ® ^K*. (l) 

where a =_L, || refers to a transversely and longitudinally polarized K* , respectively. 
(An analogous result applies to the decay into a pseudoscalar meson, but we specify our 
notation to a vector meson for simplicity.) In this equation represent universal heavy- 
to-light form factors 0, ^ and $ light-cone-distribution amplitudes. The factors Ca 
and Ta are calculable in renormalization-group improved perturbation theory. Previous 
treatments of exclusive radiative decays correspond to evaluating (|l|) to leading logarith- 
mic accuracy (up to a weak annihilation contribution as we discuss below). The main 
result of this paper is to show that systematic improvement is possible and to extend 
the accuracy to the next order. In Section ^ we present the result of the calculation. In 
Section ^ we first discuss the inputs to our numerical result, such as Wilson coefficients 
and meson parameters. We then analyse the stability of the result with respect to the 
input parameters and a variation of the renormalization scale. The phenomenological 
analysis is given in the subsequent two sections. Section ^ is devoted to S — > K*-^] Sec- 
tion 1^ to -B — > K*i^£~ . We focus in particular on corrections to the forward-backward 
asymmetry, which is independent on the form factors to first approximation. In the 
final Section ^ we discuss power corrections related to photon structure on a qualitative 
level and summarize our conclusions. 



2 Non-factorizable corrections 

We generally neglect doubly Cabibbo-suppressed contributions to the decay amplitude. 
The weak effective Hamiltonian is then given by 
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(2) 



i=l 



and we use the operator basis introduced in |^ for the operators Oj, z = 1, . . . , 6. (More 
precisely, due to the normalization of (0), Oi = 4Pj with Pi defined in 0.) We define 



Op 



9,10 



2n 



(4) 



where ttem = glml^'^ is the fine-structure constant and rhh{^) the b quark mass in the 
MS scheme. We shall later trade the MS mass for the pole mass rrib and the PS mass 
mb^PS. To next-to- leading order the pole and MS masses are related by 



rhbin) = mb 1 + 



31n- 



+ Oial 



(5) 
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(a) (b) 

Figure 1: Leading contributions to {YK*\Hcq\B). The circled cross marks 
the possible insertions of the virtual photon line. 

with as = <ys{f^)- The sign convention for Oyg corresponds to a negative Cj^s and 
+igsT^, +igemGf for the ordinary quark-gauge-boson vertex (e/ = —1 for the lepton 
fields). We will present our result in terms of "barred" coefficients Ci (for i = 1, . . . ,6), 
defined as certain linear combinations of the Ci as described in Appendix |^. The linear 
combinations are chosen such that the Ci coincide at leading logarithmic order with the 
Wilson coefficients in the standard basis ||^. 

As for form factors and non-leptonic two-body decays there exist two distinct classes 
of non-factorizable effects. (By "non-factorizable" we mean all those corrections that are 
not contained in the definition of the QCD form factors for heavy-to-light transitions. 
For example, the familiar leading-order diagrams shown in (a) and (b) of Figure |l| are 
factorizable.) The first class involves diagrams in which the spectator quark in the B 
meson participates in the hard scattering. This effect occurs at leading order in an 
expansion in the strong coupling constant only through a weak annihilation diagram 
[Figure |l|c] . The relevant diagrams at next-to-leading order are shown as (a) and (b) in 
Figure ^ below and in Figure ||. They contribute at order to the functions Ta in 
(H). Diagrams of this form have already been considered (for = 0) in §]. However, 
bound state model wave-functions (rather than light-cone distribution amplitudes) were 
used and no attempt was made to systematically expand the hard scattering amplitude 
in l/rrib. As a consequence, the result of ||^ for i? ^ K*'-f depends on an infrared cut-off. 
This difficulty is resolved in the present factorization approach. The second class contains 
all diagrams shown in the second row of Figure ^ below. Here the spectator quark is 
connected to the hard process represented by the diagram only through soft interactions. 
The result is therefore proportional to the form factor C,a and the hard-scattering part 
gives an cts-correction to the functions Ca in (|I]). 

In this section we present the results of the calculation of these diagrams. Some of the 
results needed for diagrams of the second class can be extracted from work on inclusive 
radiative decays |10[ and we have made use of these results as indicated below. The 



conventions for the form factors and light-cone distribution amplitudes for B mesons 
and light mesons are those of M. 



2.1 Notation and leading-order result 

Since the matrix elements of the semi-leptonic operators Og^io can be expressed through 
B — i> K* form factors, non-factorizable corrections contribute to the decay amplitude 
only through the production of a virtual photon, which then decays into the lepton pair. 
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We therefore introduce 



-I %{q^) [{Ml - m\,) e*^ - {e* ■ q) {p'' + p"' 

,2 



Ml - ml. 



(6) 



where H^e denotes the weak effective Hamiltonian. We define the overall quark mass 
factor as the pole mass here. The matrix element decomposition is defined such that the 
leading order contribution from the electromagnetic dipole operator Oj reads Ti{q^) = 
C-jTi{q^) + . . ., where Ti{q^) denote the tensor form factors. Including also the four- 
quark operators (but neglecting for the moment annihilation contributions), the leading 



logarithmic expressions are 

T,{q^)=CfUq^)+Y{q^) 

CfT,{q') + Y{q') 

CfUq') + Y{q^) 



%{q') 



2mh{MB + rriK* 

,2 



q 



2mi)[MB — rriK 
Mb - rriK 



2mh 



Mq 



viq% 

A,iq'), 

2, Mb + rriK* 



2mb 



A{q' 



(7) 
(8) 

(9) 



with Cf = Cr- C3/3 - 4C4/9 - 2OC5/3 - SOCg/g = Cr- {AC3 - ^5)79 - (4C4 - Ci)jZ 
and 



= m,) (3Ci + + 3C3 + + 3C5 + 

- i h[s, m,) (4 (Cs + C4) + 3^5 + a) - ^ Ks, 0) (Cs + 3^4 



2 



(10) 



The function 



h{s,mq) = - 
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^ 2 
In 



^ -9(2 + ^) 



k - II < 



arctan 



z > 1 



1 + z in 

in — z < 1 

sTz 2 - 



(11) 

is related to the basic fermion loop. (Here z is defined as 4m^/s.) Yi^s) is given in the 
NDR scheme with anticommuting 75 and with respect to the operator basis of 0. Since 
Cg is basis-dependent starting from next-to-leading logarithmic order, the terms not 
proportional to h{s,mg) differ from those given in [0. The contributions from the four- 
quark operators Oi-q are usually combined with the coefficient Cg into an "effective" 
(basis- and scheme-independent) Wilson coefficient Cg^{q'^) = Cg + Y^q"^). 
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The results of this paper are restricted to the kinematic region in which the energy 
of the final state meson scales with the heavy quark mass in the heavy quark limit. 
In practice we identify this with the region below the charm pair production threshold 
2 < 4m2 ^ 7GeV\ The various form factors appearing in are then related by 

symmetries ||. Adopting the notation of [||, simplify to 



q 



^T^{q' 
2E 



Mb 

Mb 
2E 



Cf 5i + 



q^ 



2mbMB 



2mh 



(12) 
(13) 
(14) 



where E = (M^ — q"^) /{2Mb) refers to the energy of the final state meson and ^_Ljj refer 
to the form factors in the heavy quark and high energy limit. The factors 6i are defined 
such that 6i = 1 + 0{as)- The a^-corrections have been computed in |^ and will be 
incorporated into the next-to-leading order results later on. The appearance of only 
two independent structures is a consequence of the chiral weak interactions and helicity 
conservation, and hence holds also after including next-to-leading order corrections 
T2| . We therefore present our results in terms of the invariant amplitudes 71 ||(g^), which 
refer to the decay into a transversely and longitudinally polarized vector meson (virtual 
photon), respectively. At next-to-leading order we represent these quantities in the form 

JSr Mb ^ J 00 Jo 



where Cp = 4/3, A^'^ = 3, 



-a,± 



= 1, 



[u,uj) 



= mK*/E, and Ta^±{u,u) is expanded as 



rpiO) I N I 

7;\±(u,o;) + 



u, u). 



(16) 



denotes the usual K* decay constant fx*- fK*,± refers to the (scale-dependent) 
transverse decay constant defined by the matrix element of the tensor current. The 
leading-order coefficient C^^-* follows by comparison with Eqs. (IT! 



and (|Tj) setting 5i = 1. 
To complete the leading-order result we have to compute the weak annihilation am- 
plitude of Figure |l|c, which has no analogue in the inclusive decay and generates the 
hard-scattering term Tj^°]_(m,u;) in (15). To compute this term we perform the projec- 



tion of the amplitude on the B meson and K* meson distribution amplitude as explained 
in . The four diagrams in Figure |l]c contribute at different powers in the 1 / ruh expan- 
sion. It turns out that the leading contribution comes from the single diagram with the 
photon emitted from the spectator quark in the B meson, because this allows the quark 
propagator to be off-shell by an amount of order mfjAqcD, the off-shellness being of order 
ml for the other three diagrams. With the convention that the K* meson momentum 
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is nearly light-like in the minus light-cone direction, the amplitude for the surviving 
contribution depends only on the minus component of the spectator quark momentum. 
This is in contrast to the discussion in where, using the same convention for the 
outgoing meson, the hard-scattering amplitude depended only on the plus component of 
the spectator quark momentum. As a consequence the B meson distribution amplitude 
now appears as coefficient of in the B meson projector. Except for this 
difference the calculation follows the rules of [Q. The result reads 

t|°V {u,uj) = Ti°l {u,uj)= T|[°| {u,uj) = (17) 

T| '_{u,uj) = -eg — : C3 + 3C4 . (18) 

Mbuj — — le rrib 

At leading order weak annihilation occurs only through penguin operators with small 
Wilson coefficients. For B^ decay, there is an additional Cabibbo-suppressed contri- 
bution involving current-current operators, which can be obtained by multiplying the 
previous equation with the factor 

Vus^ub Ci + 3C2 ^^g^ 



whose modulus is about 0.7. Both contributions are numerically small relative to the 
dominant leading-order terms. Note that since ^|(q'^) does not contribute to the decay 
amplitude in the limit — 0, there is no weak annihilation contribution to -B — K*'y 
at leading order in as and at leading order in l/rrib. 

We emphasize that our discussion of weak annihilation refers only to the leading 
contribution in the heavy quark limit. It is a novel aspect of 5 — > K*l^l~ (with the 
K* longitudinally polarized) that weak annihilation does not vanish in this limit, if 

~ mbAqcD- On the other hand, the sum of transverse and longitudinal contributions 
is dominated by the transverse decay amplitude at small g^, and it is probable that the 
power-suppressed annihilation contributions to the transverse amplitude 7j^(g^), which 
we neglected, are numerically more important. For = annihilation contributions of 



this sort have been studied in |13|. The interplay of the various terms merits consideration 
for b —>■ d transitions, in which weak annihilation through current-current operators is 
neither CKM- nor a^-suppressed. Since our main analysis concerns b ^ s decays, we 
defer this issue to a later discussion. 



2.2 Next-to-leading order — spectator scattering 



The hard scattering functions T^^l. in ([T5| ) contain a factorizable term from expressing the 
full QCD form factors in terms of ^q, related to the a^-correction to the 6i in Eqs. ([T^ ) 
(ll) above. We write 



(1) 



rp{i) I rp{ni) 



^ uE 



Cf 



The factorizable correction reads 
2Mf 



B 



2mbMB 



Y{q' 



2M| 



(20) 
(21) 
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Figure 2: Non-factorizable contributions to (7*^*|ifefr|-B)- The circled cross 
marks the possible insertions of the virtual photon line. Diagrams that follow 
from (c) and (e) by symmetry are not shown. Upper line: hard spectator scat- 
tering. Lower line: diagrams involving a, B ^ K* form factor (the spectator 
quark line is not drawn for these diagrams). 



Tf^_{u,uj) = Tl'i{u,uj) = Q 



(f) 



(22) 



The non-factorizable correction is obtained by computing matrix elements of four-quark 
operators and the chromomagnetic dipole operator represented by diagrams (a) and (b) 
in Figure |^. The projection on the meson distribution amplitudes is straightforward. In 
the result we keep only the leading term in the heavy quark limit, expanding the ampli- 
tude in powers of the spectator quark momentum whenever this is permitted by power 
counting. In practice this means keeping all terms that have one power of the spectator 
quark momentum in the denominator. Such terms arise either from the gluon propagator 
that connects to the spectator quark line or from the spectator quark propagator, when 
the photon is emitted from the spectator quark line. We then find: 



.(nf) 



(m, uj 



u 



+ 



B 



eut±{u, rric) (C2 + C4, - Cq 



+ Cd t±(u, nik) (C-g + ^4 - a - Amb/MB C^) + t^{u, 0) ^3 



M, Ul) 





M. 



B 



rrib 

+ 6^^11(^,0)^3 
u, u) = 



eJw {u, rric) {C2 + C4 - Cq) + t\\ {u, rrib) (C3 + C4 - Cq) 



(23) 
(24) 

(25) 



Mbuj 



Mbuj - g2 



rics 



u + uq^/M 



+ 



6M, 



B 



(h{uMl + uq^, {C2 + ^4 + C^) + h{uMl + uq\ m,) (C3 + 0^ + 0^ 
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+ h{uMl + uq\ 0) {C; + 3^4 + SCe) - — (C3 - - 15C, 



(26) 



Here Cg^^^ = + - C4/6 + 2OC5 - lOCg/S = Cg + (4^3 - ^5)73, e„ = 2/3, = -1/3 
and eq is the electric charge of the spectator quark in the B meson. h{s, mg) has been 
defined above. The functions ta{u, nig) arise from the two diagrams of Figure §3 in which 
the photon attaches to the internal quark loop. They are given by 



t±{u,mg) 



t\\{u,mg) 



2M, 



B 



uE 

2Mb 
uE 



hijng) + 
h{mg) + 



(5o(mM| + uq\ mg) - Bo{q\ nig)) (27) 
{BQ{uMl + uq\mg)- Bo{q\mg)), (28) 



where Bq and Ji are defined as 

Bo{s,mg) = - 



1 arctan 



1 



lArrig 



2ml 



and 



Li(x) 



± 



In- 



1 

X - 



m„ 



1/2 



^^\mi}J s — 1 
+ Li(x_) - Li(y+) 



-^1(1/- 



(29) 
(30) 



uM% + ug2 



y± 



2 U 



2\l/2 



ln(l 



Lio 



X — 1 



(31) 



(32) 



The correct imaginary parts are obtained by interpreting m^ as m^—it. Closer inspection 
shows that contrary to appearance none of the hard-scattering functions is more singular 
than I/m as M — >■ 1. It follows that the convolution integrals with the kaon light-cone 
distribution are convergent at the endpoints. 

The limit q^ ^ Q {E ^ Mb/ 2) of the transverse amplitude is relevant to the decay 
B K*'~f. The corresponding limiting function reads 



q)\q^=0 



U \ 



2ml 
uMl 



Li{x+) + I/i(x_; 



(33) 



In the same limit the longitudinal amplitude develops a logarithmic singularity, which 
is of no consequence, because the longitudinal contribution to the B — >• K*l^l~ decay 
rate is suppressed by a power of q^ relative to the transverse contribution in this limit. 
It does, however, imply that the longitudinal amplitude itself is sensitive to distances of 
order 1/y/q^ and not perturbatively calculable unless q^ ~ mfeAgcD- (This long-distance 
sensitivity appears already at leading order in (ffl) through the light-quark contributions 
toy(g2).) 
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2.3 Next-to-leading order — form factor correction 



The next-to- leading order coefficients C^^ in (15) contain a factorizable term from ex- 
pressing the full QCD form factors in terms of ^a, related to the a^-correction to the 5i 
in Eqs. (|l|), We write C^^^ = CP + C(°^). The factorizable correction reads 



C 



(f) 



Cf [4\n 



(f) 



mt 



4-L 



(34) 
(35) 



with 



L = 



1 1 
— q 



In 1 



mt 



(36) 



The brackets multiplying Cj^ include the term 31n(m^//i^) — 4 from expressing the MS 
quark mass in the definition of the operator in terms of the b quark pole mass accord- 
ing to (^. The non-factorizable correction is obtained by computing matrix elements of 
four-quark operators and the chromomagnetic dipole operator represented by diagrams 
(c) through (e) in Figure ||. 

The matrix elements of four-quark operators require the calculation of two-loop di- 
agrams with several different mass scales. The result for the current-current operators 
Oi^2 is presented in as a double expansion in /ml and mdmi,. Since we are only 
interested in small g^, this result is adequate for our purposes. (Note that only the result 
corresponding to Figure la-e of |]TD[ is needed for our calculation.) The 2-loop matrix 
elements of penguin operators have not yet been computed and will hence be neglected. 
Due to the small Wilson coefficients of the penguin operators, this should be a very good 
approximation. The matrix element of the chromomagnetic dipole operator [Figure 0c] 
is also given in in expanded form. The exact result is given in Appendix |B[ All this 
combined, we obtain 



(nf) 



CfC 



— C2F2 



(7) 



CfF^ 



2mbMB 



r pC^) -L r^eff p(7) 



+ 



2-^2 

Mb 
2mh 



1 



(9) 



+ CfF^^ 



(37) 



(38) 



The quantities Fi'^2^^ 



.(7,9) 



can be extracted from ||TO|-Q The quantities F^ ' ' are given m 
or can also be extracted from in expanded form. In expressing the 



Appendix 

result in terms of the coefficients Ci_2; "we have made use of fP + F^'>/6 = 0. We also 
substituted Cg by C^^, taking into account a subset of penguin contributions. 



^We thank M. Walker for providing us with additional data points which cover the range 0.25 < 
mc/fTT-b < 0.35 needed for the subsequent numerical analysis. 
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(d) (e) (f) 

Figure 3: Vertex corrections to the weak annihilation diagram in Figure |I|c. 
Only the photon coupling to the spectator quark line contributes at leading 
order in l/m^ as indicated by the circled cross. 

2.4 Weak annihilation 

In addition to the two sets of next-to-leading order corrections discussed up to now, there 
exist "vertex corrections" to the weak annihilation graph, as shown in Figure ^ The 
distinction between weak annihilation and the previous spectator scattering diagrams is 
not clear-cut. The diagrams in Figure 03 with the photon attached to any of the external 
lines can also pass as a correction to the leading-order weak annihilation graph. In fact, 
the scale-dependent logarithm in h{q'^,mq) in ( ]26|) cancels part of the scale-dependence 
of the leading order annihilation amplitude, (|l^). 

For 6 — > s transitions current-current operators with large Wilson coefficients cannot 
be contracted to give diagrams of the type shown in Figure |^. The total contribution 
is therefore suppressed by three factors: the strong coupling constant; the small Wilson 
coefficients of QCD penguin operators; and the suppression of the longitudinal amplitude 
in the decay rate for small g^. We shall see subsequently that the leading-order weak 
annihilation effect is already small. Hence, although the radiative corrections to this 
leading-order term are of some conceptual interest in connection with renormalization of 
the B meson distribution amplitude, we shall not include the tiny radiative corrections 
represented in Figure ^ into our numerical analysis. 



2.5 Summary 

We now summarize our main result. The matrix element {'y*{q, fi)K*{p', e*)\Hefi\B{p)) is 
not scheme- and scale-independent in the standard operator formalism with an on-shell 
basis, unless the photon is on-shell. Therefore ||(g^) are not physical quantities. From 
the expressions for the decay rates given below, it will be clear that the following three 
quantities are independent of the conventions chosen to renormalize the weak effective 
Hamiltonian: 
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C,4,^) . + ^ . + Y(e) + Cf + , ^ ^ (40) 

The quantity C9j|(g^) depends on the charge of the B meson through e^, but this is 
suppressed in our notation. The elhpses denote the a<j-corrections calculated above and 
defined through (|l^). The explicit expressions for the quantities defined in (|15D are given 
in (P^-(P5D and (^D|)-(^B]). We note that C7, Cg^x and Cgji depend on the conventions 



adopted in to define the soft form factors and ^y. In addition C^ depends also on 
the h quark mass renormalization scheme (but irtbC-j does not). We discuss this further 
after (|60|) below. 

We verified that our results are scale-independent to the required order. More pre- 
cisely, we have 

^lj-^[C,,C,,^{q^),C,,\\{q^)}=0{ala,C,.,). (42) 

The uncancelled terms at order as are proportional only to the small Wilson coefficients 
of QCD penguin operators. They are related to the unknown two-loop matrix elements 
of penguin operators, and to the incomplete evaluation of weak annihilation effects in the 
case of C9j|(g^). Numerically, however, the missing contributions are negligible. When 
the strong interaction corrections to (p9D-(p|) are included, the next-to-leading logarith- 



mic expression for Cf^ must be used. However, because Cg ~ \yi{Mw / jj) ~ l/a^ at 
leading order, the coefficient Cg is needed to next-to-next-to-leading logarithmic order. 
The relevant initial conditions for the renormalization group evolution can be taken from 



1^ and we have incorporated them into our analysis. (The necessary renormalization 
group formulae are summarized in Appendix y.) However, the three- loop anomalous di- 
mension matrix needed for the evolution is presently not known. Except for this (and the 
numerically insignificant terms discussed above), our result for B —>■ K*i^i~ is complete 
at next-to-next-to-leading logarithmic order. 



3 Numerical analysis of C7, Cg Cg || 
3.1 Specification of input parameters 

Wilson coefficients. We begin our discussion of input parameters with the Wilson coeffi- 
cient Cg, which we need to next-to-next-to-leading logarithmic order. In our analysis we 
use the complete next-to-next-to-leading logarithmic expression derived in Appendix |C|, 
in which we set the unknown 3-loop anomalous dimensions k^^^ and 7^^^ to zero. To 
estimate the uncertainty this introduces we note that the neglected terms read 

SC, = - i^^'^ D,{^^, M^) + «:(-^) D_i(^, M^) V ^^''^^'f ^^^i. vA , (43) 

2Po V 4/?o + 7r-7r A 
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Table 1: Wilson coefficients at the scale /i = 4.6 GeV in leading-logarithmic 
(LL) and next-to-leading-logarithmic order (NLL). Input parameters are 
Ag^ = 0.220 GeV, mt{mt) = 167 GeV, Mw = 80.4 GeV, and sin^Ow = 0.23. 3- 
loop running is used for as- We also give Cg^io at NNLL obtained as described 
in the text. 







C2 






C5 




LL 
NLL 


-0.257 
-0.151 


1.112 
1.059 


0.012 
0.012 


-0.026 
-0.034 


0.008 
0.010 


-0.033 
-0.040 




Cf 


^8 


Cg 


Cio 




,^NNLL 
•-^10 


LL 
NLL 


-0.314 
-0.308 


-0.149 
-0.169 


2.007 
4.154 




-4.261 


4.214 


-4.312 



where all notation is defined in Appendix |C| and the subscript '2' refers to the com- 
ponent of the six-component vector in brackets. Under reasonable assumptions on the 
pattern of the unknown anomalous dimension matrix (vector) we observe that the largest 
contribution to the right hand side comes from k^2^. A rough estimate is therefore 



(5C„ ^ -10-^ kP 



•-2 



0.74 



(44) 



Allowing 1^2^'' I < 100, we conclude that at scales of order of the b quark mass Cg is 
known in the Standard Model to an accuracy of about ±0.1. (If k^2^ is much smaller 
than the upper limit we assume, other terms may dominate SCg, but the uncertainty 
estimate should remain valid.) This is to be compared to an error of ±0.05 due to 
the error on the top quark MS mass, mt{rht) = (167 ± 5) GeV. The value of Cg at the 
scale /i = 4.6 GeV is given in Table ^ together with the remaining Wilson coefficients. 
The electroweak input parameters that go into these numbers are summarized along 
with other input parameters in Table below. The strong coupling is always evolved 
according to the 3-loop formula. 

Quark masses. Up to now we have assumed that is the pole mass. It is usually 
assumed that the pole mass of the b quark is less well known than the MS mass or the 
potential-subtracted (PS) mass [l^. (The PS mass replaces the pole mass in quantities in 
which the b quark is nearly on-shell, but has some large infrared contributions removed.) 
We therefore replace mb by the PS mass through the relation 

4« 

mb = nib^ps {fif) + —^fif (45) 
ovr 

and use mfe ps(2 GeV) = (4.6±0.1) GeV |T6[ as an input parameter. An exception to this 
replacement is applied to the function Y{s) in which we keep the b quark pole mass - 
computed according to ( ^Bj ) - in the small contributions from b quark loops. Apart from 
reinterpreting as mb ps(2GeV) the only consequence of introducing the PS mass is 
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Table 2: Summary of input parameters and estimated uncertainties. 





80.4 GeV 




180 ± 30 MeV 


rht{rht) 


167 ± 5 GeV 




(3 ± 1) GeV-i 


\Vts\ 


0.041 ±0.003 


/i^^±(lGeV) 


185 ± 10 MeV 




1/137 


fK*,\\ 


225 ± 30 MeV 


^^QCD 


220 ± 40 MeV 




0.2 ±0.1 


m6,Ps(2GeV) 


4.6 ±0.1 GeV 


a2(i?*)±,|i 


0.05 ±0.1 


mc 


1.4 ± 0.2 GeV 


MBi\\{S))/{2mK*) 


0.47 ±0.09 








0.35 ±0.07 



an additional term 4/ij/mb that appears in the round brackets multiplying Cf^ in (|3 
We use the charm quark pole mass, rric = (1.4 ± 0.2) GeV. 

Meson parameters. Many of the relevant meson parameters are not directly known from 
experiment. QCD sum rule calculations substitute for this lack of information and we 
use as our source of information. A summary of all input parameters together with 
their assumed errors is given in Table We now discuss some of the input parameters 
in more detail. 

By convention |Q the soft form factors (including the one for a pseudoscalar kaon) 
at zero momentum transfer are related to the full form factors by 

WO)=ff(0). UO)-jf^V^-m. ^«,(0) = <-(0). (46) 

to all orders in perturbation theory. For Aq*{0) we have taken the QCD sum rule 
result from [|17| and this gives the soft form factor ^||(0) listed in Table ^ The choice 



^±(0) = 0.35 ± 0.07, also given in that Table, requires comment since (|46|) and V^* (0) 
from |T^ would give ^±(0) = 0.39 ± 0.06. The motivation for our choice of input is 
related to the fact that the QCD sum rule prediction for Tf*{0)/V^*{0) deviates from 
the behaviour expected in the heavy quark limit ||^. An alternative way of computing 
^±{0) uses ^±{0) ^ 0.93Tf *(0) (with Tf*(0) evaluated at the scale ^ = rub). Taking 
Tf" (0) from fl^ gives ^±(0) ~ 0.35. Using this smaller input has the advantage that 
our result for the B —>■ K*'~f decay rate is automatically consistent with an alternative 
representation of this decay rate in terms of the full QCD form factor T/^*(0). (This will 
be discussed in more detail later.) Furthermore, we shall see that the decay rate predicted 
for B —>■ K*'j at next-to-leading order favours form factors smaller than those estimated 
with QCD sum rules. We may then take the point of view that ^±{0) is determined by 
r[B — > K*'y) and hence also T/^*(0) and V^*{0) through the relations of This logic 
of course implies that we take seriously the heavy quark limit prediction for the decay 
rate and form factor relations in the heavy quark limit. The energy dependence of the 
form factors is assumed to be 
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as predicted by power counting in the heavy quark hmit. 

The hard-spectator scattering (and annihilation) contribution depend on the hght- 
cone wave functions of the 5-meson and the hght meson in the final state which are 
characterized in terms of decay constants, Gegenbauer coefficients etc. The leading 
contribution in the heavy quark mass expansion involves the leading-twist distribution 
amplitudes of light mesons only. We truncate the expansion of these functions into 
Gegenbauer-polynomials Cp^^^ at second order, i.e. (a =±, ||) 

^K'M = 6w(l + a,{K*)a Cf'^'\2u - 1) + a2{K*)a C^'^'\2u - 1)} . (48) 



The values for the Gegenbauer coefficients are taken from |T^, but for simplicity we 
neglect the small differences between ai{K*)±^ and ai{K*)\\. We also enlarge the er- 
ror, see Table ^ In the same Table we also give numerical values for the light meson 
decay constants |T^. (The Gegenbauer coefficients and the decay constant for a trans- 
versely polarized vector meson, f±, are scale-dependent and assumed to be evaluated 
at the scale 1 GeV. The variation of Gegenbauer moments has a negligible effect on 
our result compared to the overall parameter uncertainties. We therefore neglect the 
scale-dependence of the Gegenbauer moments in the numerical analysis, but we evolve 
f± according to /±(/u) = /±(/Uo) («s(/u)/as(;Uo))^^^^.) The renormalization scale in the 
hard-scattering and annihilation terms is chosen lower than in the form factor contribu- 
tions in order to reflect the fact that the typical virtualities in the hard scattering term 
are of order mbAqco rather than m^. If /ii (assumed to be of order m^) is the scale 
in the form factor term, we choose (/iiA/^)^/^ with A/j = 0.5 GeV in the hard-scattering 
term. This convention applies to all scale-dependent quantities in the hard-scattering 
term including and the Wilson coefficients. 

The two B meson light-cone distribution amplitudes enter only through the two 
"moments" 

A-=r,.£s±M, (49) 

JO UJ 

xtA<f)=r^ ^f;;^' . (50) 

Jo oj — /Mb — le 

The moment of is identical to the moment that appears also in non-leptonic B 

decays [^, the decay B [|l^ and in the factorization of form factors [Q. The 



appearance of ^b- at leading order in the heavy quark expansion is a new aspect of the 
decay B ^ K*£+i-. 

The light-cone distribution amplitudes $b,±(cc;) appearing in (^9]), ( pO[ ) are not well- 
constrained at present and provide a major source of uncertainty in our calculation. 
Some general properties have been derived in the literature [§, The equations of 
motions for the light quark in the 5-meson lead to 

$B,-(C^) = / —<^bA^/v) ^ <^B,+ iuj) = -UJ<l>'jj_{uj) (51) 

JO rj 

We conclude from this that A^^_,_ = $^,-(0). Furthermore A^^_(g^) must diverge loga- 
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Figure 4: The absolute value of the ratio of i?-meson moments | A^^_ (g^) | / A^^_,_ 
as a function of q^. The S-meson distribution amplitudes are taken as in (|53D 
with ojq^ = (3±1) GeV-^ 

rithmically for — > 0, provided that A^^^ ^ 0: 

limA^:_(g^) = A^:+ /' ^ - oo (52) 

q2^0 Jo T] 

On the other hand, for = O(MsAqcd) the moment A^^_(g^) is finite and of or- 
der 1/Aqcd- The equations of motion for the heavy quark relate the first moments {u) b,± 
to the mass difference Ahqet = Mb — rrib, leading to {oj)b,+ = '2{uj)b,- = 4Ahqet/3. For 



Xb,+ the upper bound 4Ahqet/3 has been derived in [jT8[. The simple model functions 



= 4 e'"^"" ' ^^bA^) = — (53) 



with luq = 2Ahqet/3, that have been proposed in []Ty[, are consistent with these general 
constraints, leading to \b,+ = 2Ahqet/3, and 

>^bUq^) = Ei(gVMBCUo) + tn (54) 

where Ei(z) is the exponential integral function. Numerically, for Ahqet — 500 MeV, 
one obtains A^^^_ ~ 3 GeV~^ with an estimated error of about 1 GeV~^. For the same 
parameter values the absolute value of A^^_(g^), normalized to A^^^ is plotted in Figure |[ 
There exist alternative proposals for 5-meson light-cone wave functions, for instance 
from the Bauer-Stech-Wirbel model pO| or variants of it. Since in these models the 
distinction between $B,+ (t^) and ^b,-{^^) is not made, we refrain from presenting a 
thorough comparison of different models and stick to (^4]) to evaluate the moment of 
$B,_(a;). This means that we neglect a systematic uncertainty related to the shape of 
but this uncertainty is irrelevant numerically, because the moment in question 
appears only in the small annihilation contribution and the small correction to Cg j|(g^). 

All input values from the meson sector together with their estimated uncertainties are 
summarized in Table |^. We note that apart from the renormalization scale uncertainty 
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Figure 5: Renormalization scale-dependence of jCyp at leading (LO) and next- 
to-leading order (NLO). The curve "NLOi" shows the NLO result without the 
spectator scattering correction. 



and the error in our knowledge of a^, the most important uncertain parameters can be 
collected into a single factor 

fBfK*,a .ggX 

N,MbXb,+ UO) ^ ' 

that determines the relative magnitude of the hard-scattering versus the form factor 
term. Adding all errors in quadrature, this factor is uncertain by about ±50%, where 
the largest error is currently from A^^^. 

3.2 Exclusive effective "Wilson" coefficients 

Having specified our numerical input, we now discuss the three effective "Wilson" coef- 
ficients C7, C9,_L(g^) and C9j|(g^). 

We begin with the quantity to which the decay rate of 5 — > K*'-) is proportional. 
In Figure |^ we show the renormalization scale dependence of this quantity at leading 
and at next-to-leading order. We also show a curve that corresponds to setting the hard- 
scattering term to zero, i.e. to taking into account only the correction C^^-* in (pIS]). The 



reason for considering this term separately is that it should cancel the sizeable leading- 
order scale dependence, while the hard scattering correction is a physically different 
effect that appears first at next-to-leading order. Figure § shows that this is indeed 
correct. The hard scattering correction reintroduces a mild scale-dependence. The most 
important effect is however a large enhancement of jCyp at next-to-leading order. At the 
scale rrib = 4.6 GeV we find 

ICtInlo/ICtIlo ^ 1-78, (56) 

which corresponds to a sizeable, but not unreasonable 33% correction on the amplitude 
level. The form factor and hard-scattering correction contribute about equally to this 
enhancement. More precisely, the [non-]factorizable part of C^'^ (defined in (p!5D) is a 
—8% [+24%] correction to the real part of the amplitude, the [non-]factorizable part 
of t\^\ is a +11% [+5%] correction (all numbers at /x = mb). The error on jCyp is 
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estimated by combining in quadrature the error from all input parameters as specified 
in the previous section and from allowing the renormalization scale to vary from mb/2 
to 2mb. (The scale in the hard scattering term is accordingly lower, see the discussion 
above.). The result is 

I^tInLO = ^■^'^^-0.026^ (57) 

where the largest errors are from scale dependence (±0.014), A;^^^^ (±0.015), ^^(0) 
(±0.009), and Aqcd (±0.010). The value given refers to using the PS mass mb_ps(2 GeV), 
since only the product mfe^ps(2GeV) \Cr\ is convention- independent. The NLO result for 
must be compared to ICyIlq = 0.098 with a 30% error from scale dependence alone. 
At this place it is also appropriate to note that our result is based on the heavy quark 
expansion, and therefore there exist corrections of order Aqcd / '"^6 at the amplitude level. 
At present we have no means of quantifying these corrections in a systematic way. 

The g^-dependent coefficients Cg^±{q'^), C9j|(g^) defined in (^), (^1]) are shown in 
Figure The left panels display the reduction of renormalization scale dependence in 
going from leading order to next-to-leading order (including the intermediate approxi- 
mation with the hard scattering term set to zero). The right panel gives the complete 
NLO result with all input parameter uncertainties and scale dependence included in the 
error estimate. 

The characteristic features of the result can be understood from the discussion of C7 
above. Since Cg > and C^^ < 0, the contribution to Cg^ail'^) from the virtual photon 
matrix element, Ta{q'^), is always negative (a =±, ||). In the case of a transverse virtual 
photon, this contribution is enhanced by a factor M^/g^ owing to the real photon pole 
and hence dominates Cg^ j_ at small q^. This causes Re (Cg^x) to change its sign and results 
in the characteristic shape of |Cg_x|. Due to the large contribution from the photon pole 
the next-to-leading order contribution is again substantial and the value of momentum 
transfer q^ at which Re(Cg^x) changes its sign is increased. Since q^ determines the 
location of the forward-backward asymmetry zero, this important correction will be 
discussed in some detail in Section The photon pole is absent in the longitudinal 
coefficient Cgjj, so the NLO correction to it is much smaller. The form factor and 
hard scattering correction (both small) nearly compensate each other leaving no net 
modification of the leading order coefficient. Contrary to the transverse coefficient, the 
longitudinal one depends on the charge of the spectator quark in the B meson, i.e. on 
the charge of the decaying B meson. This effect is already present at leading order via 
the annihilation contribution. At NLO it amounts to a few percent of Cg_ y , the precise 
magnitude being g^-dependent. We return to a discussion of this isospin-breaking effect 
in the context of the B — > K*i~^i~ decay spectrum. 

The peculiar discontinuity in the error band around = 6 GeV^ in two of the 
right panels of Figure || is unphysical and related to the fact that in computing the 
error we allow the charm quark mass to be as small as 1.2 GeV. The discontinuity is 
a consequence of the cc threshold in the charm loop diagrams and reminds us that the 
validity of our calculation is not only limited by the requirement that the K* momentum 
is large (which puts an upper bound on g^), but also by the lack of a model-independent 
treatment of the resonant charmonium contributions. The physical "threshold" begins 
at Mj^^ ^ 9.6 GeV^, but the perturbative approximation is expected to fail earlier. 
Model-dependent studies of the charmonium contributions suggest that the perturbative 
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Figure 6: Momentum-transfer dependence of jCg^j,!, jCg, \\ \ and Re {Cg^±). The 
left panels show the LO, NLO and NLO result without spectator scattering 
("NLOi") for fixed input parameters. The width of the bands demonstrates 
the renormalization scale dependence estimated by variation from mb/2 to 
2mb- The right panel shows the next-to-leading order result including renor- 
malization scale and input parameter uncertainties (all added in quadrature), 
and for comparison the leading order result (dashed). The solid line in the 
center of the band refers to the NLO result with default parameters. jCg ||| is 
given for B~ mesons. 
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approximation should be valid up to g2 ~ (6 - 7) GeV^ 13 . (This supposes that 
we discard a charm mass as small as 1.2 GeV in estimating effects related to the charm 
threshold.) 



4 The decay B K*j 

The decay rate for B — > K*'y in the heavy quark limit is given by 



^TT^ ^ \ Mi An 



with mf, PS = m;, ps(2 GeV) the PS mass. We included the kaon mass squared in a phase 
space and kinematic correction, but we generally neglect such terms in the dynamical 
quantity 7^(0) = ^±(0)C7. By helicity conservation the kaon is transversely polarized 
and only the transverse soft form factor ^_l(0) appears in the decay rate. In numerical 
form, the branching fraction is 



Br(.. K^,) ^ m ■ (^) (If )^ ^ (7.9-;^) ■ 10^ 



(59) 



using \C'j\^ = 0.175l|]:l]i^ and m6,ps = (4.6 ± 0.1) GeV, ^±(0) = 0.35 ± 0.07, \Vt 



ts 



0.041 ± 0.003, = 1/137, as in Table ^. (r^ denotes the B meson lifetime.) For 
comparison we note that the leading order prediction is 4.5 ■ 10~^. The large difference 



2 

reflects the enhancement of \Ci\ at next-to-leading order as discussed in Section 3.2 



The branching ratio predicted for the default parameter set is nearly twice as large as 



the current experimental averages |23, 24, 25 



Br(5° ^ is:*%)cxp = (4.54 ± 0.37) ■ 10-^ (60) 
Bi{B- i?*"7)cxp = (3.81 ± 0.68) • 10"^ (61) 

Note that we predict no difference of the neutral and charged B meson decay rates at 
leading order in the heavy quark expansion, so that the main effect of isospin breaking 
on branching fractions is the different lifetimes of the charged and neutral mesons. 

Before speculating on the origin of this discrepancy, we note that there are alternative 
ways of representing the result of our calculation. Since the decay amplitude is propor- 
tional to a single form factor, T/^*(0), there is no need to introduce the soft form factor 
^±(0), and we can express the decay rate directly in terms of the full QCD form factor. 
(The need to introduce the soft form factors arises when we consider the K*i~^i~ final 
state, which involves many form factors that all relate to the same soft form factors.) 
Furthermore the dominant dependence on the b quark mass arises through the factor 
rrib in the definition of the operator Ot, so it also more natural to keep the MS mass 
m}){rh}y) as a prefactor rather than the PS mass. This alternative representation amounts 
to a redefinition of Cy in which all factorizable corrections vanish at the scale rhf,, if we 
understand T/^*(0) as renormalized at the scale rhb. We have computed the coefficient 
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C'y also in this alternative scheme (indicated by the prime) and find that (|58D is modified 
as follows: 



\C,f = 0.175l°:°l \C',\' = 0.165t°-°^^ 



(62) 



-0.017 

The error is smaller in the primed scheme, because the hard-scattering correction is 
much reduced once the factorizable correction is reabsorbed into the full QCD form 
factor. This implies less sensitivity to the uncertain parameters that normalize this 
correction. The two values given in ( |62D are consistent with each other if one relates 
Tf *{0) = 1.08 ^x(O) and takes into account that mfe,ps(2 GeV) = 4.6 GeV corresponds 
to Tfibirhb) = 4.4 GeV. 

Several facts may explain the difference between the predicted and observed decay 
rate. The first and most interesting possibility is a modification of the standard model 
at short distances that would result in a smaller value of IC7I. However, the modification 
needed to explain the observed decay rate is far too large not to be ruled out by the 
agreement between experiment and theory for the inclusive decay B — >■ Xg'y. It appears 
equally implausible that new interactions would modify only the spectator scattering 
and hence not show up in the inclusive rate. We therefore conclude that the explanation 
must be sought in our understanding of QCD effects. Most of the NLO enhancement 
is related to the non-factorizable form factor type correction which appears in a nearly 
identical form in the inclusive decay rate as well. Hence this enhancement cannot be 
simply dismissed without putting the agreement for the inclusive decay into question. 
A second and less interesting possibility is therefore to invoke a large l/rrih correction 
that would render our calculation unreliable. Given the smallness of the non-factorizable 
hard scattering correction it is not obvious how such a large dynamical enhancement of 
l/ruh terms could be explained. This is in particular so as large l/rrib corrections known 



to exist for non-leptonic decays 0, ^ such as "chirally enhanced" corrections and large 
weak annihilation contributions are absent for decays into vector mesons. 

We therefore consider seriously the possibility that the form factors at = are 
substantially different from what they are assumed to be in the QCD sum rule approach 
or in quark models. This possibility is entertained in Figure |^, where we show how 
the experimental data constrains the soft form factor ^_l(0) and the normalization of 
the hard scattering term. (We take all parameters other than A^^^_ in the normalizing 
factor constant and let A^^_|_ be a representative. That is, a variation of A^^^_ should 
be considered as a variation of a combination of parameters of the type (0). All other 
parameter uncertainties are approximately accounted for by adding a 20% uncertainty 
to the measured branching fraction which we combine with the experimental error in 
quadrature.) The result of this fit is that 

T/^*(0)|^=^, = 0.27 ±0.04 [a(0) = 0.24 ±0.06], (63) 

where the error may be loosely interpreted as a one standard deviation error. This 
determination yields a smaller form factor than a similar determination in |]12|, where 
the NLO correction is not included and the bottom quark mass is handled differently. 

Our results can be extended in a straightforward way to the decay B p^y. This 
decay is particularly interesting in searches for extensions of the standard model, because 
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Figure 7: Fit to C-l(O) and \b,+- The solid curves give the la, 2a and 3a 
regions that follow from the observed B — > K*'y decay rate. The light solid 
curves mark the ranges assumed for ^±(0) and Xb,+ in Table ^ 



of the suppression oi h ^ d transitions in the standard model and the simultaneous 
chirality suppression. Except for trivial adjustments the most important difference to the 
K*'-) final state is that the different flavour structure allows weak annihilation to proceed 
through the current- current operator with large Wilson coefficient C*2. This annihilation 
contribution is power-suppressed in 1/mft, but the suppression is compensated by the 
large Wilson coefficient and the occurrence of annihilation at tree level. The annihilation 



contribution is calculable and has to be taken into account in a realistic analysis p3| , p^ . 
This will be done elsewhere p8|, and we restrict ourselves to 6 ^ s transitions in 
this paper. 



5 The decay B K*i+£- 



For the decay B K*i~^i we obtain the double differential decay spectrum (summed 
over final state polarisations, and lepton masses neglected) 



dq^dcosO 1287r3 



M%\{q\m 



2 ^2 ^3 ^ Clom 

47r 



K* 



i + cos^^)|^e±(gy(|C9,x(g^)r + c: 



Ml 



+ cos^ 



'mi 

, rriK* 



2\2 



cos^^^e±(g')'Re(C9,^(g2))Cio 



(64) 



where 



A(g ,mj^,) 



.2 \2 



1 - 



Ml, 



Ml 



+ 



K* 



Ml) Ml 



1/2 



(65) 
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and is the invariant mass of the lepton pair. The angle 9 refers to the angle between 
the positively charged lepton and the B meson in the center-of-mass frame of the lepton 
pair. We have kept terms of order m\* in an overall factor related to phase space 
and kinematics but emphasize that the expression in brackets neglects such terms. The 
first two terms with angular dependence (1 ± cos^ 6) correspond to the production of 
transversely and longitudinally polarized kaons, respectively. The third term generates 
a forward-backward asymmetry with respect to the plane perpendicular to the lepton 
momentum in the center-of-mass frame of the lepton pair. The factor Ay multiplying the 
Wilson coefficient Cio in (0) arises from the factorizable corrections to the form-factor 
A2 and is given by (cf. Eq. (66) in 0]) 



"^^ ^ 47r ^ ' An NMb{E /mK*)i\\{q^) Jo u ^ ^ ^ 

with L defined in (p^) and E = (M^ — q^)/{2 Mb)- (This factor could be eliminated by 
choosing a factorization convention different from (^61) for the soft form factor and by 
redefining Cgji accordingly.) 

Without going into much detail, we remark that the result for the decay into pseu- 
doscalar mesons, B —>■ Ki~^i~, is easily obtained from the corresponding expressions 
for the decay into longitudinally polarized vector mesons. The angular distribution is 
predicted to be proportional to (1 — cos^ 9) and for the lepton invariant mass spectrum 
we obtain 

dT _ ^^l^*:^*^lV|A(g^m^)^f^Yep(gy (|C9,p(g^)P + C,^o) (67) 



where the quantities analogous to those in (0), (40) are now defined as 



and 



{j%q,f,)Kip')\H,,\Bip)) = ^^y:svj-^^^[q'{v'+v')-m-<)ci' 

(69) 

The non-factorizable contributions to Tp are the same as those to T\\ up to an overall sign 
and the replacement 'E.\\fK,\\^K*,\\{u) — >• fK^xiu) in (p!5|). The factorizable contributions 
follow from the ratio of the tensor form factor friq^) and the soft form factor ^p{q^) as 
calculated in The net result for Cg^p is then 

where the replacements Sy/^, ||$_ft'*,||(M) — fx^Kiu), ^\\{q'^) — ^ — ^p(q'^) are understood 
to be performed on the right-hand side. 
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Figure 8: Differential decay rate dBT{B~ —>■ K*~i^i~)/dq'^ at next-to-leading 
order (solid center line) and leading order (dashed). The band reffects all the- 
oretical uncertainties from parameters and scale dependence combined, with 
most of the uncertainty due to the form factors ||(0). 



5.1 Lepton invariant mass spectrum 

The lepton invariant mass spectrum dBr/dq"^, obtained by integrating (16^) over cosO, 
is shown in Figure (We use the B meson lifetime tb- = 1.65ps.) As the decay rate 
is dominated by the contribution of longitudinally polarized K* mesons and the contri- 
bution from Oio except for small g^, the impact of the next-to-leading order correction 
is small for > 2 GeV^. (This can be directly inferred from the small correction to 
ICgjil as seen in Figure ^.) The apparently rather large uncertainty of our prediction is 
mainly due the form factors with their current large uncertainty and to a lesser extent 
due to \Vts\ and the top quark mass. It may be hoped that in the longer term future 
the form factors could be known with much greater confidence. The uncertainties can 
then be reduced by a factor of three in which case the limiting factor of our prediction 
is probably l/rrib corrections, all other parameter and scale variations being negligible 
in comparison. 

Since \Cg^ y | depends on the charge of the decaying B meson, there is a small amount 
of isospin breaking. While this effect is g^-dependent for jCg yl, the combination of coeffi- 
cients that enters the lepton invariant mass spectrum turns out to be almost independent 
on q^. The effect is, however, very small, 

dTjB- K*~i+£~)/dq^ - dT{B^ K*H+£-)/dq^ 
~ dT{B- K*-i+i-) /dq^ + dT{BO K*H+i-)/dq^ ^ ^ 

For the differential branching fractions the main isospin breaking effect arises from the 
lifetime difference of the charged and neutral B mesons. 

The impact of the NLO correction on the lepton invariant mass spectrum to the decay 
in a pseudoscalar meson follows a qualitatively similar pattern as for the vector meson 
final state, in particular as the decay rate involves only a quantity closely related to \Cg^ \\ \ 
(see ([70|)). We therefore do not discuss this decay in further detail here, but note that 
the lepton invariant mass spectrum develops a logarithmic singularity for small g^. This 
is due to the long-distance sensitivity mentioned above, which is now dominant because 
the photon pole is absent. The invariant mass spectrum is therefore non-perturbative 
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Figure 9: Forward-backward asymmetry dAp-B{B~ K*" l'^ / dq^ at next- 
to-leading order (solid center line) and leading order (dashed). The band re- 
flects all theoretical uncertainties from parameters and scale dependence com- 
bined. 



for ~ ^QCD5 but perturbative for q^ ~ m^AQCD- Furthermore, the non-perturbative 
contribution is formally power-suppressed when the lepton invariant mass spectrum is 
integrated from to some q^ of order m^AQCD- 



5.2 Forward-backward asymmetry 

The QCD factorization approach proposed here leads to an almost model-independent 
theoretical prediction for the forward-backward asymmetry It has been noted in 



|31| that the location of the forward-backward asymmetry zero is nearly independent of 
particular form factor models. An explanation of this fact was given in where it 
has been noted that the form factor ratios on which the asymmetry zero depends are 
predicted free of hadronic uncertainties in the combined heavy quark and large energy 
limit. In |Q the effect of the (factorizable) radiative corrections to the form factor ratios 
has been studied and has been found to shift the position of the asymmetry zero about 
5% towards larger values. We are now in the position to discuss the effect of both, 
factorizable and non-factorizable radiative corrections to next-to-leading order in the 
strong coupling constant on the location of the asymmetry-zero, and hence to complete 
our earlier analysis. 

We define the forward-backward (FB) asymmetry (normalized to the differential de- 
cay rate dV{B- K*-(i+(i-)/dq^) by 

dA-pB 



dq" 




(72) 



Our result for the FB asymmetry is shown in Figure ^ to LO and NLO accuracy. From 
(|6^) it is obvious that dA-pB/dq^ oc Re (C9^j_(g^)), and therefore the FB asymmetry van- 
ishes if Re (C9_^(gQ)) = 0. At leading order this translates into the relation 

C, + Re{Y{ql)) = -^-^^Cf , (73) 

% 
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which, as aheady mentioned, is free of hadronic uncertainties. This effect can also be 
seen in Figure ^ where the uncertainty for the prediction of the FB asymmetry gets 
smaller around q^. Using our default input parameters we obtain ql = 3.4l!;o;5 GeV^ 
at LO where the error is by far dominated by the scale dependence. The numerical 
effect of NLO corrections amounts to a substantial reduction of the FB asymmetry for 
intermediate lepton invariant mass (g^ = 2 — 5 GeV^) and a significant shift of the 
location of the asymmetry zero to 



ql = 4.39t°:i GeV^ 



(74) 



The largest single uncertainty (about ±0.25 GeV^) continues to be scale dependence 
(varying within mh/2 < fi < 2mb as usual), but also the variation of m^, Aqcd, and 
\b,+ gives an effect. 



Our LO result for differs somewhat from the value gg = 2.8811q_28 quoted in 



The reason for this difference is a combination of three effects: (i) slightly different 
values for the SM Wilson coefficients, in particular a larger value for Cg in [^, (ii) a 
different treatment of the 6-quark mass, (iii) the inclusion of (higher-order factorizable) 
form factor corrections in 



since the full QCD form factors estimated from light- 
cone QCD sum rules fl^ are used. The first two issues are resolved within the present 
next-to-leading order treatment but with respect to (iii) we note that the form factor 
ratios entering Cg^± are exactly those where there exists a discrepancy between the QCD 
sum rule result and the heavy quark limit 0|. The origin of this discrepancy has not 
been clarified, but it is conceivable that it is related to some 1/Mb corrections, which are 
implicitly included in the sum rule estimate. In order to estimate the impact of using full 
QCD form factors on the location of the asymmetry zero, we reabsorb the factorizable 
as corrections into the physical form factors and the MS mass. The FB asymmetry then 
reads (with the tensor form factors and the MS mass rhf, all renormalized at the scale /i) 



dA 



FB 



dq^ 



1287r3 



dV/dq^ 

C, + Y{q') + 



MlKq' 



m 



K* 



Up 
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8f_ 

Ml 



CroA^{q')V{q' 



xRe 



2^ , OtsCp Q{nl,9)i^_2 



Air 




+ ^ (Ms + rriK* 



q^ y ■ V{q^) 



(MB-mK^) 1 



MIJ A^{q^ 



du^K' ±(m)t[°2(m 



(75) 



47r Nc Mb 

The non-factorizable contributions from C^^^^ in (pT]) are now divided into two parts 
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coming with the operator structure of and Og, respectively, 



CFCf''\q') = - [CFr'iq') + CrFr{q')\ • (77) 

These are just the terms that also enter the inclusive decay, while the last term in ([75|) 
arises from the hard spectator-scattering and thus only appears in the exclusive decay. 
Using the relations T2 = 2ETi/Mb and Ai = 2EMbV/ (Ms+mft-.)^, which are presumed 
to hold to all orders in as, and to leading order in l/mj, 0, ([75| ) can be rearranged 
to contain only V and Ti/V. If furthermore the latter ratio is expanded in as and if 
some terms of order m\*/M'^ are neglected, we return to the next-to-leading order result 
that we use by default. On the basis of ([75|) and the form factor estimates from [l7], we 
now obtain ql = 3.25 GeV^ at LO and gg = 3.94 GeV^ at NLO. (To obtain this result 
we use rrib^ps as input and treat perturbatively the difference with rhb-) The NLO value 
is slightly more than one "standard deviation" away from our default result in (|7^). 
The difference between the two values could be viewed as an estimate for an additional 
systematic error, but the main conclusion is that using the soft form factors or the full 
QCD form factors from light-cone QCD sum rules does not cause a large ambiguity in 
our result for ^q. We therefore assume 

q^ = (4.2 ± 0.6) GeV^ (78) 

to be a conservative estimate for the asymmetry zero. 

The Wilson coefficients may also receive contributions from new particles in theories 
beyond the standard model. (Of course extensions of the standard model may introduce a 
larger set of operators as well.) The function Cg.j, depends dominantly on Cg and Cj^ (see 
(^)). The experimental measurement of the decay rate for the inclusive decay B — * Xs'j 
already constrains \Cj^\ to be close to its standard model value. The forward-backward 
asymmetry will establish the sign of C^^ and assuming this to take its standard model 
value, a measurement of the FB asymmetry zero in the decay B K*i^i~ provides a 
way to determine Cg. It is therefore instructive to consider our result as a prediction for 
Cg for a given value of gg that will be measured in future experiments. In Figure |10| we 
show our NLO prediction for Cg at the scale /i = 4.6 GeV as a function of q^ together 
with the standard model prediction (see Section For comparison we also show the 
LO result and the modification arising from using ( [75| ) instead of (0). We remark that 
Cg is scheme- dependent and so the plot refers to the particular renormalization scheme 
for Og adopted in ^j. 

Figure |10| is our main phenomenological result. It illustrates that the FB asymmetry 
zero can be used to test the standard model and to search for extensions, since the 
theoretical uncertainties relevant for the determination of Cg are under control. In this 
respect the systematic inclusion of NLO corrections, which is provided in this work, turns 
out to be essential, because it reduces the renormalization scale and scheme dependences 
to a large extent and corrects the leading-order result by a large amount. We emphasize 
that in contrast to the absolute rates for B K*'y and B — > K*l^l~ the dependence 
on the B — * K* form factors is sub-leading for the FB asymmetry zero. Thus the 



cM'\q') + [F^^q^) + (g^)) + Cf (g^) 



(76) 
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Figure 10: The Wilson coefficient C^^imi^ deduced from a measurement of 
the asymmetry zero go at next-to-leading order including all parameter un- 
certainties (light band), but not including the slightly enlarged error in ([78|). 
For comparison we also show (without uncertainties) the leading order result 



(dashed) and the result at next-to-leading order using (^) and the full QCD 
form factors from [|T^ (long-dashed). The horizontal band is the standard 
model value C9(4.6 GeV) = 4.21 ± 0.12. 



current difficulty to account for the B K*^ branching ratio (which may suggest 
a smaller value of the form factor Ti(0) than usually estimated from QCD sum rules 
as discussed above), does not necessarily cause a problem for our prediction of Cg vs. 
^Q. We also note that the resulting values of are sufficiently below the threshold for 
charmonium resonances {J/ip-, so that a long-distance contribution to the function 
y (g^) can safely be neglected. We therefore conclude that if the asymmetry zero is found 
at some Qq < 6 GeV^, the short-distance coefficient 6*9(4.6 GeV) can be determined with 
a theoretical error of ±10%. If its value turns out to be different from the standard 
model value 4.21 ± 0.12 by an amount A and if we assume that all Wilson coefficients 
other than Cg^io remain unmodified at the scale Mw, then (|97[) below implies a new 
contribution of magnitude A to the coefficient function Cg at the scale Mw- 



6 Concluding discussion 

Building on our previous work on heavy-to-light form factors |Q and on work on non- 
leptonic decays [@, Hi, we demonstrated that also the radiative B decays B V'j and 
B Vi^i~ can be computed in the heavy-quark limit. This allows us to solve the prob- 
lem of non-fact or izable strong interaction corrections (i.e. those corrections not related 
to form factors), which has so far prevented a systematic discussion of exclusive radiative 
decays, as compared to their inclusive analogues. The approach presented here does not 
circumvent the need to know the heavy-to-light form factors, but we may hope that 
progress in lattice QCD will give us these form factors at smaller g^ and more reliably 
than at present in the longer term future. 

In the present paper we have concentrated on the b s transition. We noted that the 
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next-to-leading order correction yields an 80% enhancement of the B K*'y decay rate, 
a large part of which is related to the NLO correction that appears also in the inclusive 
decay |10|- The enhancement is so large that the predicted decay rate disagrees with 



the observed decay rate unless the form factors at = are much smaller than what 
they are believed to be, or unless our theoretical prediction is made unreliable by a large 
correction to the heavy quark limit. Perhaps the most interesting outcome of our analysis 
is that while the NLO correction to the lepton invariant mass spectrum in B ^ K*£~^£~ is 
very small, there is a large correction to the predicted location of the forward-backward 
asymmetry zero. The particular interest in this quantity derives from the fact that 
all dependence of form factors arises first at next-to-leading order. This implies that 
(given the Wilson coefficient C7) the Wilson coefficient Cg can be determined from a 
measurement of the location of the asymmetry zero with little theoretical uncertainty. 
We find that this conclusion holds true in particular after including the NLO correction, 
which however, shifts the predicted location of the zero by 1 GeV^, or about 30%, with 
a residual uncertainty estimated to be about (0.4 — 0.6) GeV^. This implies that an 
accurate measurement of the asymmetry zero determines the Wilson coefficient Cg at 
the scale rrib with an error of ±10%. 

The method developed here can also be applied to b —>■ d transitions. These are 
particularly interesting in the search for modifications of the standard model. The better 
control of standard model effects after accounting for non-factorizable corrections should 
help increasing the sensitivity to such non-standard effects. In particular we note that 
the approach presented here allows us to compute isospin breaking effects, which may 
turn out to be a particularly nice signal of non-standard physics. A detailed discussion 
of b ^ d transitions requires, however, a more careful study of weak annihilation effects, 
and we plan to return to non-standard physics in the context of such a study pU[ . 

A systematic, but rather difficult to quantify, limitation of the factorization approach 
is the poor control of l/rub suppressed effects. These arise from various sources and we 
conclude by mentioning one of them. We may ask to what extent a non-leptonic decay 
B K*p in which the p meson is assumed to convert to a photon through vector meson 
dominance may affect our calculation of the direct decay into K*'y*. More accurately, 
the question is to what extent the hadronic structure of the photon matters at small q^. 
It is not difficult to see that the indirect contribution is suppressed by one power of in 
the heavy quark limit. For the sake of the argument we return to the vector dominance 
picture and note that there is a one-to-one correspondence between diagrams that appear 
in the factorization approach to non-leptonic decays p, |^ and the diagrams computed 
here. For example, the naive factorization contribution to the former case corresponds 
to our leading order weak annihilation term; the spectator scattering contribution to 
non-leptonic decays corresponds to the diagram in Figure |^ with the photon attached 
to the internal quark loop and so on. Counting powers of rrib for all the quantities that 
appear in these amplitudes we find that the direct amplitude is always larger by a factor 
of nif,. 
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Note added 



While completing this work, the article appeared, in which the next-to-leading order 



hard-scattering correction io B fq is calculated. The corresponding quantity should 
be obtained from the function T^^^^ in (|23|) in the limit 0. The small QCD penguin 
contributions have been neglected in , so that only the terms proportional to C2 and 
have to be retained for the comparison. The charm quark mass in t±{u,mc) is set 
to in presumably because an endpoint divergence is found for rric 7^ 0. We note 



that our result does not exhibit this endpoint divergence (which would invalidate the 
factorization approach), but setting rric = for the sake of comparison, we obtain 



u Znih u 

We then differ from (4.12) of by a sign in the first term and a factor of —4 in the 
second term. We also note that includes a factorizable form factor correction from 
^ while retaining the full QCD form factor ^1(0). This implies that the factorizable 
correction is included twice. 

A Operator bases 

We use the renormalization conventions of [0, but present our result in terms of the 
following linear combinations of Wilson coefficients: 

C, = ici, 

C2 = C2 — - Ci, 

D 

C'3 = C3-^C4 + 16C5-^C6, 

6 

C*4 = - C4 + 8 Cq, 

C5 = C'a — - C4 + 4 C5 — - Cq, 

O 6 

C, = ^C, + 2C,. (79) 

These definitions hold to all orders in perturbation theory. The "barred" coefficients are 
related to those defined in by 
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B Unexpanded form of 



(7,9) 



Using the conventions of |IU| we obtain for the 1-loop matrix element of Os (Figure ^ 
together with a set of symmetric diagrams): 



32 , fi 8 s , . 8 4 11 - 16s + 8s2 
— In In s m ; — ; 

9 rrib 9 1-s 9 9 (1 - s)^ 

4 1 



.(9) 



9 (l-s)3 

16 1 , . 

m s + 



(9s - Ss^ + 2s^) Bo{s) - (4 + 2s) Co(s) 
8 5-2s 8 4-s 



[(l + s)fio(S)-2Co(s) 



9 1 - s ■ 9 (1 - s)2 9 (1 - s)3 
where s = q'^/ml, Bo{s) = BQ{q'^,ml) is given in (p9D , and the integral 

1 



Cn(s) 



dx ■ 



In- 



X 



10 x(l — s) + l 1 — x(l — x)s 
can be expressed in terms of dilogarithms. 



(82) 
(83) 



(84) 



C NNLL formulae for Cg and Ci ... ^ 

Because the anomalous dimension of Og begins at order q;°, the Wilson coefficient Cg is 
needed to next-to-next-to-leading logarithmic (NNLL) accuracy. This requires also the 
coefficient of the four-quark operators to this accuracy. Although not all of the necessary 
inputs are currently known, we present here the solution of the renormalization group 
equations to NNLL order. 

We first restrict ourselves to the sector of four-quark operators. Define as(/i) = 
Q!s(/i.)/(47r) and (3{as) = Pof^l + ■ ■ ■ with /5o = 11 — 2nf/3. We further define the matrix 
U (/i, Mw) such that 

C{fx) = Uifx,Mw)CiMw), (85) 
where C(/i) is the vector of Wilson coefficients. U{^, Mw) satisfies 



d 

dlnfi 



Uifx,Mw)=l''ifJ^)Uifx,Mw) 



(86) 
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and the anomalous dimension matrix is expanded as 
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Let V be the matrix that diagonahzes 7^°) , so that 



and define 



w / 



(0) 



diag 



V 



-7fV(2/3o)' 



V 



(87) 



(88) 



(89) 



diag 



which solves (p6| ) to leading order in as(/i). Then the NNLL expression for the evolution 
matrix U (/x, Mw) reads 



where 



The matrices if*^^) have the entries 



(1) 



(0) A 



(1) 



ij 



2P^ 2/30 + 7^-7 



JO) ' 



(90) 
(91) 

(92) 



(2) 



. (0) /32 , ^ 2/3o + 7f ^ - 7f 

"^l^i AQ2^ 2^ .n , (0) (0) 



(1) rrCi) _ ^//(i)^ 



ij "jk 



G 



(2) 



4/3o + 7f^-7 



.(0) 



(93) 



and we have defined 



(94) 

The same equations can also be used to solve the renormalization group equations when 
the electromagnetic and chromomagnetic dipole operators are included. Then 7 is an 
8x8 matrix, and one needs to keep in mind that the 6x2 block in 7*^") that mixes 
the dipole and four-quark operators is obtained from n + 2 loop diagrams, while all 
other entries come from n + 1 loop diagrams. In our numerical analysis we expand 
U{fi, Mw) C{Mw) in powers of as{fi) and as{Mw) using ( pOD and the expansion of the 
initial condition C{Mw) in powers of as{Mw) and keep only terms to the required order 
in both expansion parameters. 

The renormalization group equation for Cg is particularly simple and can be obtained 
directly by integrating the coefficients of four-quark operators. The equation reads 



d 

dlnfj, 



(95) 
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where C (/i) is the vector of Wilson coefficients of four-quark operators and 



(96) 



is the 1x6 matrix that describes the mixing into Og. The calculation of k^"-* involves 
n + 2 loop diagrams. The solution of is given by 



CM = C,{Mw) + Mw) C{Mw), 



with the 1x6 matrix 



ly (/i, Mw) = -- / da, f/(/i, Mw) 

2 Jas(Mw) p[as) 



(97) 



(98) 



with the evolution matrix f/(yU, My/) from the four-quark sector as defined above. The 
solution to NNLL accuracy is obtained by inserting ?7(/i, My/) to this accuracy. We 
introduce the 6x6 matrices 



D^{^i,Mw) = V 



1 



n-7fV(2A)) 



V-\ (99) 



diag 



in terms of which the solution is given by 



W{^^,M, 



w) 



2/3oasiMw) 



D_i{fx,Mw) 



1 

Wo 
a. (Ml 



K 



2/3o 



- ^ + J^A Doifi, Mw) - /^^-^^ D_,{f,, M^) J^'^ 



\Po PoJ V I3q J 

+ J(2) j Di(/x, Mvk) - - |i + k^-^) Z^)^ Do(/^, M^^) 



ft 



(-1) Z^_i(/x,Mvk) 



(100) 



In our numerical analysis we expand (|97D in powers of as{My/) and keep only terms to 
the required order in the expansion parameter. To NNLL we thus need Cq{Mw) to order 
Qg and since the ffist term in the expansion of this quantity is of order a^, we recover 
from (|100|) the well-known result that the LL solution for Cg is induced only by mixing 
and is of order l/a^. 
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